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Abstract 

In this paper we consider a stepping-stone model on a circle with circular Brownian 
migration. We first point out a connection between Arratia flow and the marginal 
distribution of this model. We then give a new representation for the stepping-stone 
model using Arratia flow and circular coalescing Brownian motion. Such a represen- 
tation enables us to carry out some explicit computation. In particular, we find the 
Laplace transform for the time when there is only a single type left across the circle. 
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1. Introduction 

Stepping-stone model is a mathematical model for population genetics. A discrete- 
site stepping-stone model describes the simultaneous evolution of interacting popula- 
tions over a collection of finite or countable colonies. There are mutation, selection 
and resampling within each colony, and there is migration among different colonies. 
See |Kim53j and |Shi88j for some early work. 

Continuous-site stepping-stone model was first introduced in |Eva97j . It is a process 
takeing values from the space 

where E denotes the continuous site space, IK denotes the type space, and M%(K) 
denotes the space of all probability measures on IK. Intuitively, such a map \x simulta- 
neously represents the relative frequencies of different types in populations at various 
sites. More precisely, for e 6 E and B C K, /j,(e)(B) represents the "proportion of 
the population at the site e processing types from the set B" . The "moments" of 
the continuous-site stepping-stone model are specified using the so called migration 
processes taking values in E. 

In this paper we only consider a stepping-stone model with site space T, a circle of 
circumference 1, with type space K = [0, 1], and with Brownian migration on T. We 
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call it a stepping-stone model with circular Brownian migration (in short, a SSCBM) 
and write it as X throughout the paper. 

The distribution of SSCBM is uniquely determined by a family of coalescing Brow- 
nian motions on T. But we have to go through more notations before we could present 
the explicit formula. 

Given a positive integer n, let V n denote the set of partitions of N n := {1, . . . , n}. 
That is, an element n of V n is a collection ir = {A\, . . . ,Ah} of disjoint subsets of 
N n such that (Ji^-i = ^n- The sets A±, . . . Ah are the blocks of the partition it. The 
integer h is called the length of ir and is denoted by |7r|. Equivalently, we can think of 
V n as the set of equivalence relations on N n and write i ~ n j if i and j belong to the 
same block of ir £ V n - 

Given ir £ V n , let 

ai := min^j, 1 < i < \ir\. 

{oi} is the collection of minimal elements for ir. 

By a circular (instantaneously) coalescing Brownian motion we mean a collection 
of Brownian motions on T such that any two of them will move together as soon as 
they first meet. Given a circular coalescing Brownian motion (Z\, . . . , Z n ) starting 
at e = (ei, . . . , e n ). For t > 0, let ir e (t) be a ^-valued random partition such that 
i ~7re(f) j iff Zi(t) = Zj(t). Then vr e (t) is the random partition induced by (Zi(t)). 
Write 

T e (t) :={ ai (t):l<i<\7T e (t)\} 

for the collections of minimal elements for vr e (t). 

SSCBM is then a H-valued Hunt process X with its transition semigroup specified 
as following. Given fj, G 3 and n > 0, for any fi £ C(T), C K, % = 1, . . . , n, 

(1.1) 



n „ 

TT / deif(ei)Xt(ei)(Ki 
i=i 



where denotes the probability law of X when its initial value is [i. See Theorem 



4.1 in DEF + 00 for a result on a general continuous-site stepping-stone model. 

In |DEF + 00] a particle representation for X was given using the Poisson random 
measure on Z?x[0, oo[x!K and a "look down" scheme similar to that in D K96j . It leads 
to better insight into the model. In the same spirit we are going to propose another 
representation for X in this paper. 

It was shown in |Eva97j that X degenerates, i.e. for any t > 0, for almost all e£l, 
Xt(e) becomes a point mass on some k € K. A stronger version of this clustering 
behavior was later shown in |DEF + 00] for site-space T and in |Zho03| for site-space 
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R. In fact, when the site-space is R there exists a random partition of R such that 
R is divided into disjoint intervals and X t (e) is a point mass on the same k £ K for 
almost all e in each interval. This suggests that we can identify X t with a function / 
on T such that /(e) = k whenever Xt(e) = 5k- In this way we can identify X with a 
step-function-valued process. 

Let H' be the space of K- valued right continuous step functions on T equipped with 
the topology inherited from Dr(T). For each /i £ S, we are going to construct a 
H'- valued process {X[,t > 0) which can be regarded as SSCBM with initial value \x 
under the above-mentioned identification. To this end, we first point out an interesting 
connection between Arratia flow and SSCBM in Section 2. This connection allows us 
to specify the entrance law of X' using the pre-image of Arratia flow. Then we give 
an explicit construction of X' using Arratia flow and circular coalescing Brownian 
motion, and we will show that X' so defined does have the right distribution under 
the above-mentioned identification. In this sense X' provides a nice version for X. 
Such a representation enables us to compute the distribution of the time when there is 
only a single type of individuals left across T in Section 3. It also allows us to obtain 
a result on the type that survives eventually. 

2. A REPRESENTATION OF STEPPING-STONE MODEL WITH CIRCULAR BROWNIAN 

MIGRATION 

We adopt some conventions for the rest of this paper. We identify T with in- 
terval [0,1). Whenever we write (ei,...,e m ) £ T m it implies that ei,...,e m have 
been already arranged in anti-clockwise order around T. Given u, v £ T, write [u,v[ 
for an interval starting at u and ending at v in anti-clockwise order. Write v — u 
for the length of the interval [u,v[. For {ki} C IK and (ei,...,e m ) £ T m , write 
YllLi [}, e m+ i := ei, for a right continuous step function on T. 

Arratia flow was first introduced in |Arr79j . Arratia flow on T describes the evo- 
lution of a stochastic system in which there is one Brownian motion starting at each 
point in T. Two Brownian motions coalesce once they meet. Formally, the Arratia 
flow can be defined as a collection {(p(s,t,x) : < s < t, x £ T} of random variables 
such that 

• the random map (s,t,x) i— > <p(s,t,x) is jointly measurable, 

• for each s and x, the map 1 1— > (j>(s,t,x), t>s,is continuous, 

• for each s and t with s < t, the map x i— > (p(s,t,x) is non-decreasing and 
right-continuous , 

• for s < t < u, cj)(t, u, •) o (p(s, t, •) = 4>(s, u, •), 

• for u > 0, (s, t, x) i— ► 4>(s + u,t + u,x) has the same distribution as cj), 
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• for (xi, . . . , x rn ) G T m the process (0(0, t, xi), . . . , 0(0, t, x m )) t >o has the same 
distribution as a circular coalescing Brownian motion starting at (x±, . . . , x m ). 

From the continuity of Brownian sample paths we see that, for each t > 0, there 
exists a positive integer valued random variable N(t) and two sequences of random 
variables (Vi(t)) G T N ® and (ZT^i)) G T N ® such that 

0(0, t,x) = Vi{t) forxG [t/i(t),r/i+i(*)[ and % = 1, . . . ,JV(t). 

In fact, we can even show that 

oo 

P[|iV(t)|] = 1 + 2^exp{-nVt} , 

n=l 

where |iV(t)| denotes the cardinality for N(t). See Corollary 9.3 in |DEF + 00] . 

For any fj.eE, given N(t), (Ui(t),i = 1, . . .,N(t)) and (Vi(t),i = 1, . . .,N(t)), let 
{/Cj,i = 1, . . . , iV(i)} be a collection of independent IK-valued random variables such 
that Ki follows the distribution fi(Vi(t)). Define 

N(t) 

(2.1) X[{e) = KiHWt), U i+1 (t)[}(e), e G T. 

i=l 

We first point out that X'(t), when identified as 

N(t) 

(2-2) ^^l^W^+i (*)[}, 

i=l 

is indeed a version of X t . 

Proposition 2.1. For any t > 0, wrf/i t/ie identification K2.ty X[ has the same dis- 
tribution as Xt under Q^. 

Proof. To determine the distribution of X[ we only need to specify joint distributions 
such as 

F{X' t (e 1 )edk 1 ,...,Xl(e n )edk n }. 

By definition (0(0, t, e\), . . . , 0(0, t, e n )) is a circular coalescing Brownian motion 
starting at (ei, . . . , e n ). Let 7r e (t) be the induced partition on N n . For any ki G K, i = 
l,...,n, given (Ui(t)) and (Vi(t)) as before, observe that 0(0, i, ej) and 0(0, t,ej) 
belong to the same interval [U r (t), U r+ i(t)[ for some r iff i ~ir e (t) 3- Also notice that 
0(0, t, e^ = 0(0, t, ej) whenever i ~ vr e( 4 ) j. Then 



(2.3) Fl^iXKe^edh} 



,i=i 



An inspection of (|2.3|) reveals that (|1.1() holds for X'(t) when it is regarded as E- valued. 
So X t and X 4 ' have the same distribution. 
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□ 

By Proposition 12. II we may and will suppose that Xt,t > 0, is H'-valued in the rest 
of the paper. 

We can read off some properties for Xt, t > 0, immediately from Proposition 12.11 
First, with probability one Xt (as a function of e) can only take finitely many different 
values from K. Moreover, if /i(e) is a diffuse measure for almost all e G T, then with 
probability one Xt takes different values over different intervals on T, i.e. X(e) = k 
for e G [ei,e2[ whenever X{e\) = k = X(e2). Such properties are also discussed in 
Section 10 of |DEF+00| . 

Conditioning on X s = Y1T= l ^l{[ M «i u i+i [}> P-ip shows that, given t > s, X t can 
only take values from {fcj}. Moreover, for any {kj,j = l,...,re} C {k{\ and any 
(zj) E T n , by ()1.1[) we can further show that 



f}{Xt(z j ) = fy 



(2.4) 



3=1 



X s = kil{[ui,u i+ i[} 



1=1 



f]{X s (Z 1 (t-s)) = k' J } 



where (Zj) is a circular coalescing Brownian motion starting at (zj). 

To describe the evolution of X over time we need a Lemma on duality between two 
circular coalescing Brownian motions. 

Fix y = (yi,...,y m ) € T m and z = ( Zl ,...,z n ) G T n . Let (y x ,...,y m ) be an 
m-dimensional circular coalescing Brownian motion starting at y. Let (Z\, . . . , Z n ) 
be an n-dimensional circular coalescing Brownian motion starting at z. Put 

I£(t,z) := l{Yi(t) G [zj,z j+1 [} 

and 

I£(t,y) := l{y,G[^(t),Z J+1 (t)[} 
for 1 < i < m and 1 < j < n. Recall that z n+ \ := z\ and Z n+ i := Z\. 

Lemma 2.2. The two (m x n)- dimensional arrays (I^~(t,y)) and (I^(t,z)) have the 
same distribution. 

Proof. We can prove Lemma 12.21 in the same way as Theorem 2.1 in |Zhoj . i.e. we 
first show that the corresponding duality holds for circular coalescing random walks, 
and then apply time-space scaling to obtain the desired result for circular coalescing 
Brownian motions. 

□ 
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By Lemma 12.21 we can easily derive the following side result concerning a dual 
relationship for Arratia flow. Such a result was pointed out in |Arr79j for coalescing 
Brownian flow on the real line. 

Proposition 2.3. When identified as point processes, (Ui(t)) and (Vi(t)) have the 
same distribution for any fixed t > 0. 

Proof. For any (zj) E T 2n , let (Zj) be a coalescing Brownian motion starting at (zj). 
Consider a sequence of circular coalescing Brownian motions {(<j)(0,t,xf l ))^l 1 ,m = 
1,2,.. .} such that the set {x™,i = 1, ... , m} of starting locations approaches to a 
dense set in T as m — > do. Such a sequence provides an "approximation" for the 
Arratia flow on T. Then by Lemma 12.21 



(2.5) 



n 



0(0, t,x?) \J(z 2j -i,z 2j ) 



p n 



i=l 



cr?{J(Z2j-i(t),z 2j (t)) 



On the one hand, taking limits on both sides of 1)2. 5 Jl as m — > oo, we can show that 

os(*)>n u^i-i,^-) =^J =p |n {^i-iw = z 2j (t)}^ . 

On the other hand, {Ui(t)} n (z2j-i, ^2j) = iff O^-i, 22 j) C Ui+i(t)[ for some 

i iff 4>(0,t, z 2 j-i) = Vi(t) = (j)(0,t,z 2 j) for some i. Consequently, we also have 

n j I n I 

{Ui{t)}f] \J(z 2j ^,z 2j ) = 0| = P| f|{0(O,t,z 2j _i) = 0(O,t,3jy)} I • 

Therefore, (Ui(t)) and (V^(i)) have the same avoidance function. So, the assertion 
of this Proposition holds (see Theorem 3.3 of Kal76 ). □ 

Let us go back to the stepping-stone model. We first consider a special initial value 
\x. Given 

m 

v = ^2S ki l{[ui,u i+ i[} e E, 

i=l 

write Y = (Yi) for an m-dimensional circular coalescing Brownian motion starting at 
u := (ttj) and define 



Xl = J2kiH[Yi(t),Y i+1 (t)[},t>0, 

i=l 

with the convention that l{[y,y[} := 0. 

Lemma 2.4. X' has the same distribution as X under Q" . 
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Proof. {X' t ) is clearly a Markov process from its definition. 

Given {k[, . . . , k' n } C {k i7 i = 1, . . . , m} and (vj) G T n , let Z = (Zj) be a circular 
coalescing Brownian motion starting at v := (vj). Set 



Lemma l'2 . 21 yields that 



j=l i-.k^k'j 



•aim 



i=l J [j=l L 4=1 

'b(Y(t);v)] 

[s(u;Z(t))] 

n ( m 

f]l^2k i l{[u i ,u l+1 [}(Z J (t)) = k' J 

j=l li=l 

f){X>( Zj (t))=k'^. 

It then follows from Proposition 12.11 and (|2.4j) that X' and X have both the same 
initial value and the same transition semigroup. So, they have the same distribution. 

□ 

Now we are ready to construct a representation for X with a general initial value 
(jl E 3. Given e > 0, as in (J2.H) put 

JV(e) 

(2.6) = J>il{[^(e),C/ m (£)[}. 

i=l 

Given iV(e), (Ui(e), . . . , U N ^(e)) and . . . , rejv(e))) wr ite (li) for an A r (e)-dimensional 
circular coalescing Brownian motion starting at (t/j(e)). We further define 

AT(e) 

(2.7) X t ' = '£ / K i l{[Y i (t-e),Y i+1 (t-e)[},t>e, 

i=l 

again, with the convention that l{[y, y[} := 0. Combining Proposition ^. H and Lemma 
12.41 we can easily obtain the following result. 

Theorem 2.5. Given (j, € 3 anrf e > 0, (X^,i > e) has the same distribution as 
(X t ,t> e) under 

Remark 2.6. The representation ()2.7j) suggests that SSCBM can also be thought of 
as a multi-type, nearest-neighbored voter model on T. See Chapter V in Lig85 for 
discussions on voter model. 
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Remark 2.7. A similar representation can be found for a stepping-stone model with 
Brownian migration on R. We leave the details to the readers. 

3. The first time when there is only a single type left 

In this section we are going to study properties of X using the representation given 
in Section 2. 

Treating (X t ,t > 0) as 3'- valued, put 

T := inf{* > : 3k G K,X t (e) = k,\/e £ T}. 

T is then the first time when a single type of individuals prevail all over T. It is easy 
to see from the representation (J2.7|) that 

Q^{T < oo} = 1, 

for all fj, E 3. Now we are going to find the exact distribution for T. 

We start with a preliminary result which is interesting in its own right. Let (Yi) be 
an m-dimensional circular coalescing Brownian motion starting at (t/j) € T m , m > 2. 
Let 

T m := inf{t > : Y 1 (t) = ... = Y m (t)}. 
Proposition 3.1. Given any positive integer m > 2, we have 

(3.i) p[ e -*r»] = y 

U sinh(TA) 

Proof. For z = 1, . . . ,m, write for the time when Y^+i first reaches Y{ from the 
clockwise direction. As usual, we define Y m+ \ := Y\. Since (l^+i — Yi)/y/2 is again 
a Brownian motion which starts at (yi+i — yi)/y/2 and stops whenever it reaches 
or 1/V2, Si is then the first time that the Brownian motion (Yi+i — l^)/\/2 reaches 
1/V2 before it reaches 0. We thus have 



(3.2) P[e' 



-a&i _ sinh((y i+ i -yj)\/A) 



sinh(v / A) 
See Exercise II.3.10 in |KY91j . 
Our key observation is that 

m 

P{T m < t} = (J{Si < t}, 

i=l 

and the events on the right hand side of this equation are disjoint. So, ()3.1|) follows. 

□ 



Standard argument gives the following result. 
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Corollary 3.2. For any positive integer m > 2, we have 

1 1 m 

nT m ] = 4-4 - yrf- 

8=1 

Consequently, P[T m ] attains its maximum 1/4— l/4m 2 iff all the initial values yi, . . . , y m 
are equally spaced on T. 

Remark 3.3. An explicit expression for the distribution of T m can also be found. By 
Theorem 4.1.1 in |Kni81j . we have 

HSi <t} = V2(y l+1 -yi) + -J2 sin ( V2n7r(y i+1 - yi )) exp {-nVt} . 

n=l 

Therefore, 

^ m 00 (_-|\n , . 

(3.3) P{T m <^} = v / 2 + -^^ - — - sin ^n7r(y i+ i - exp {-nVi} . 

i=l n=l 

We expect that P{T m < t} also reaches its minimum when y±, . . . ,y m are equally 
spaced on T. But we do not have a proof yet. 

Let k be the type of individuals left after time T. Then 

k = lim Xt(e), Ve E T. 

Theorem 3.4. Given \i E E suc/t i/tai //(aj) zs a diffuse probability measure for almost 
all x £T, then the Laplace transform for T has the expression 

(3.4) Q M [e -AT ]= A>Q 

sinh(VA) 

Moreover, 

(3.5) Q m {k edk}= / defi(e)(dk). 
Proof. Given m and (y^) E T m , we first observe that, by (|3.1|) . 



(3.6) lim P[e-^] = lim V " ^ = 

m ^°° m ~^°°J^l sinh(vA) sinh(VA) 



as 



max (y i+ i - y t ) -» + 

l<i<m 



Recall from the representation (|2.7j) that, for e > 0, 

N(e) 

Xi = J2 K i 1 {[Zi(t-e),Z i+1 (t-e)[}, t>e, 
i=i 
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where, given N(e), (Zi) is a circular coalescing Brownian motion starting at (C/j(e)) £ 
T N ( £ ). Put 

T(e) := inf{t > : Z 1 (t) = ... = Z N{e) (t)}. 

Notice that, given N(e), k%, . . . , /%Ye) are au different since fi(e) is diffuse for almost 
all e£T, Then e + T(e) is also the first time when X'(e) assumes a single value in IK 
for all e£T, 
Put 

A(e) := max (U i+1 (e) - U^e)). 

l<i<N(e) 

It is evident from the definition of Arratia flow and the representation (|2,7|) that 
A(e) — ► in probability and 

F{T(e) > 0} -> 1 as e + . 

In addition, Q^{T > 0} = 1. It follows from Theorem E3] and (|3~o]) that 



(3.7) 



lim ] 



e -A( £ +T( £ )). r(e) >Q 



sinh(i/A) 

Finally, by we have 

lim Q^{X t (e) E dA:} = lim F\u(Z(t))(dk)], 
where Z is a circular Brownian motion starting at e G T. (|3.5[) thus follows. 



□ 



Remark 3.5. Notice that the distribution of T does not depend on fj, as long as /i(x) 
is diffuse for almost all x £ T. 

Remark 3.6. Let X be a Brownian motion starting at < x < 1/V2. Put 

r x := inf ji > : X t = or 1/^} . 
We observe that, for the /U in Theorem 13.41 



(3.8) 



r[e -X T] = Um smh(sV2A) 
z-o+ xV2sinh(\/A) 

= hm — b — - — — - 

*-0+ p{r x . = l/y/2) 



lim P 



= 1/V2 
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Using ()3.8|) and Theorem 4.1.1 in |Kni81j we can further find an explicit expression 
for (F{T < t}. For t > 0, 

1 ( 2 °° (— l) n \ 
QP{T <t} = lim —=- y/2x + - V — sin(\/2n7rx) exp {-nVtl 

(3.9) V n=1 7 
= l + 2^(-l) n exp{-nVi}. 

n=l 

Remark 3.7. It is not hard to see from the proof for Theorem 13 . 41 that the distribution 
()3.9|) coincides with the distribution of the time when the image of Arratia flow on T 
first becomes a set of a single element, i.e. the distribution of 

r := inf{t > : (j)(0,t,x) = </>(0,t,y),Vx,y G T}. 

Again, for the n given in Theorem 13.41 for any e > 0, let interval [U^,U"[ be the 
unique interval [Ui(e), Ui + \{e)[ in (|2.6j) such that k = kf, i.e. [U' e , U"[ is the collection 
of sites at time e whose type eventually prevails. 

Proposition 3.8. For the n given in Theorem \3-4\ as e — > 0+ both (U' e ) and (V") 
converge in distribution to a uniform distribution on T. 

Proof. Clearly U" — U' e — ► in probability. Therefore, we just need to show that for 
any [a, b[c T, 

(3.10) lim ¥{U' e G [a, b[} = b- a. 

To prove ()3.1(Jj) . we first notice that, given N(e), events {Si < oo}, 1 < i < N(e), are 
all disjoint, where Si is defined as in the proof for Proposition bvrt for a coalescing 
Brownian motion starting at (J7j(e), C/j+i(e)). Consequently, by (|3.2j) 

F{U' £ G [a, &[} = P j |J {Ui{e) G [a, b[, S t < oo} 

= E< Yl (U i+1 (e)-Ui(e))l{Ui(e)e[a,b[} 

^l<i<N(e) 

= E[U n (e)-U m (e)}, 

where m := min{i : Ui(e) G [a,b[} and n : = max{i : Ui{e) G [a,b[}. Therefore, (|3.10f) 
follows readily. 

□ 



12 



Remark 3.9. If fi £ S is arbitrary, we can not find the explicit distribution for T under 
. Nevertheless, similar to the proof for Theorem 13.41 we can still show that 

oo 

Q M {T < t} > 1 + 2 ^(-l) n exp {-n\ 2 t} , t > 0. 

n=l 
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